Product of Ginibre matrices: Fuss-Catalan and Raney distributions 
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Squared singular values of a product of s square random Ginibre matrices are asymptotically 
characterized by probability distributions Ps{x), such that their moments are equal to the Fuss- 
Catalan numbers or order s. We find a representation of the Fuss-Catalan distributions Psix) in 
terms of a combination of s hypergeometric functions of the type sFs-i- The explicit formula derived 
here is exact for an arbitrary positive integer s and for s = 1 it reduces to the Marchenko-Pastur 
distribution. Using similar techniques, involving Mellin transform and the Meijer G-function, we 
find exact expressions for the Raney probability distributions, the moments of which are given 
by a two parameter generalization of the Fuss-Catalan numbers. These distributions can also be 
considered as a two parameter generalization of the Wigner semicircle law. 



INTRODUCTION 



Random matrices of various ensembles find numerous 
applications in several fields of statistical physics. In the 
general class of non-hermitian random matrices aii im- 
portant role is played by the Ginibre ensemble [1]. A 
matrix G of size N of such an ensemble consists of 
independent random complex numbers, drawn accord- 
ing to the Gaussian distribution with zero mean and a 
fixed variance [1, Q . Such matrices are used to describe 
non-unitary dynamics of chaotic systems and open quan- 
tum systems [jj. This ensemble of random matrices can 
also be used to analyze the human EEG data [S], for 
telecommunication applications based on the scattering 
of electromagnetic waves on random obstacles [6] , or in 
mathematical finances to describe correlation matrices of 
various stocks 0, [1| ■ 

The spectrum of a non-Hermitian matrix G belongs to 
the complex plane. Spectral density of random matrices 
of the suitably normalized Ginibre ensemble is described 
by the Girko circular law Q, as in the limit N ^ oo 
it covers uniformly the unit disk. The random matrix 
W = GG^ , called a Wishart matrix, is positive. Hence its 
eigenvalues Ai, i — I, . . . ,N, are real and non- negative. 
Introducing a rescaled eigenvalue, x = NX one can show 
that in the limit of the large matrix size the spectral 
density P{x) converges to the Marchenko-Pastur (MP) 
distribution (Toj . 

In general, products of random matrices are a sub- 
ject of an intensive research for many years [ll| . Recent 
studies on products of Ginibre matrices concern mul- 
tiplicative diffusion processes il2i] . correlation matrices 



used in macroeconomic time series 13| , a random matrix 



approach to quantum chromodynamics [T3] and lattice 
gauge field theories [15]. Properties of the complex spec- 
tra of products of random Ginibre matrices were recently 
analyzed in ; 16] . 

It is also interesting to study singular values of a prod- 
uct of s independent Ginibre matrices, X — Gi ■ ■ ■ Gg. 
Note that a squared singular value of the product X 
equals the corresponding eigenvalue of the Wishart like 
matrix W = XX'^. For s — 2, positive random ma- 
trices of the form W2 = GiG2(G'iG2)^ found their ap- 
plications in finances [l^. Matrices of the form Ws = 
Gi ■ ■ ■ Gs{Gi ■ ■ ■ GsY for an arbitrary s were used to de- 
scribe random quantum states associated with certain 
graphs [l3l and quantum states obtained by orthogonal 
measurements in a product of maximally entangled bases 

The corresponding asymptotic level density distribu- 
tion Ps{x) is called Fuss-Catalan distribution of order s, 
since its moments are given [19- 21] by the Fuss-Catalan 
numbers [H, [2^, (also called Fuss (sOj numbers [111)- 
Strangely enough, the Fuss-Catalan numbers generalize 
the Catalan numbers, although the work of Fufi [l^ was 
done much earlier then the contribution of Catalan 
The Catalan number can be defined as a number of dif- 
ferent bracketing of a product of n -I- 1 numbers, or the 
number of possible n folding of a map which contains 
n -\- 1 pages in a row 
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The Fuss-Catalan distribution describes asymptoti- 
cally statistics of singular values of the s-th power of 
random Ginibre matrices. This result obtained recently 
by Alexeev et al. [l^ was derived by estimating the mo- 
ments of the distribution of squared singular values of 
a power ff* of a random matrix and showing that these 
moments converge asymptotically to the Fuss-Catalan 
numbers. This is true under rather weak assumptions: all 
entries of the matrix G are independent random variables 
characterized by the zero mean, variance set to unity and 
finite fourth moment. 
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The Fuss-Catalan distribution can be considered as a 
generalization of the MP distribution, which is obtained 
for s = 1. Moreover, the distribution Ps{x) belongs to the 
class of free Meixner measures |29l|, and in terms of free 
probability theory it appears as the free multiplicative 
convolution product of s copies of the MP distribution 
[HIIoll, which is written as Ps{x) = [Pi{x)]^' . 



An explicit form in the case s — 2 was derived in 31 1 
in context of construction of generalized coherent states 
from combinatorial sequences. The spectral distribution 
of Pg {x) for a product of an arbitrary number of s random 
Ginibre matrices was recently analyzed by Burda et al. 
(3^ also in the general case of rectangular matrices. The 
distribution was expressed as a solution of a polynomial 
equation and it was conjectured that the finite size effects 
can be described by a simple multiplicative correction. 
Another recent work of Liu at al. 33] provides an integral 
representation of the distribution Ps{x) derived in the 
case of s square matrices of size iV, which is assumed 
to be large. However, these recent contributions do not 
provide an explicit form of the distribution Ps{x). 

The aim of this note is to derive exact and explicit 
formulae for the Fuss-Catalan distribution Ps{x), which 
can be represented as a combination of s hypergeometric 
functions. The derivation is presented in Sec. II, while 
some auxiliary information on special functions and the 
proof of positivity of Ps{x) are provided in Appendices 
A and B, respectively. In section III we discuss a certain 
two-parameters generalization of Fuss-Catalan numbers. 
As these numbers quantify generalized Raney sequences 
[2^ . the corresponding probability measures Wp^r{x) will 
be called Raney distributions. As special cases they in- 
clude Marchenko-Pastur distribution, Fuss-Catalan dis- 
tributions and the Wigner semicircle law. We find exact 
expressions for Raney distributions corresponding to in- 
teger parameter values in the general case and provide 
explicit formulae in the simplest cases of small values of 
the integer parameters 1 < r < p. 



II. FUSS-CATALAN DISTRIBUTIONS 

For any integer number s one can use the binomial 
symbol to define a sequence of integers denoted by 

FCs{n), 



FCs{n) 



1 



1 



S7i + n 
n 



(1) 



Here n = 0,1, . . . , while s = 1, 2, . . . , and these numbers 
are called the Fuss-Catalan numbers of order s [13] . 

We enumerate some of these sequences below for n — 
0,...,7: 

FC'iin) = 1,1,2,5,14,42,132,427,... 
FC2in) = 1,1,3,12,55,273,1428,7752,... 
FCsin) = 1,1,4,22,140,969,7084,53820,... 
FCi{n) = 1,1,5,35,285,2530,23751,231880,.... 



The above sequences are contained in the Online Ency- 
clopedia of Integer Sequences (OEIS) [131 under the la- 
bels (A000108), (A001724), (A002293) and (A002294), 
respectively. These sequences can be considered as a 
generalization of the sequence FCi{n), which consists of 
Catalan numbers, C (n) — ^^q^j (^^) . 

We are going to show that for any given s there exists 
a density distribution Ps{x), which satisfies 



x'^Ps{x)dx ^ FCs{n), n==0, 1,.. 



(2) 



where 



(3) 

In other words we are looking for a positive density Pg {x) 
which satisfies the above infinite system of equations. As 
the density turns out to be defined in a finite segment 
0, Kg], the solution of this Hausdorff moment problem 
35| associated with the Fuss-Catalan numbers is unique. 
An explicit proof of positivity of Ps{x) is provided in 
Appendix B. 

We employ the method of the inverse Mellin transform, 
which was previously used to construct explicit solutions 
of the Hausdorff moment problem [H, 113] and to derive 
explicit form of the Levy-stable distributions [1^. The 
Mellin transform of a function f{x) and its inverse 
A4^^ are defined by a pair of equations 



ria) := M[fix);a]^ f 
Jo 



x''-^f(x)dx (4) 



and 



f{x) := M-'[r{o);x] = ^. I 



c+ioo 



x-^r{a)da, (5) 



with complex a. In ^ this variable is integrated over 
a vertical line in the complex plane [s^. For discussion 
concerning the role of c, see ( [39l44l| ). Therefore the 
solution of Eq. ^ can be obtained by extending integer 
variable n to complex cr by a substitution n — !■ a — 1. The 
desired form of the distribution Ps{x) can be formally 
written as an inverse Mellin transform, 

Psix) = M-^[FCs{<j);x] . (6) 

To find such a transform we will bring the Fuss-Catalan 
numbers into a more suitable form. Representing the 
binomial symbol in ([T|) by the ratio of Euler's gamma 
function one obtains 



FCsia) 



r[(5 + i)(a-^)] 



(7) 



r[.(.-^)] r(a) 

After applying twice the Gauss-Legendre formula (|A1[) 
for multiplication of the argument of the gamma function 
one arrives at 

-1 - 



FCJa) 



1 rs + 1 



/27r 



s-3/2 



+ 



s+1/2 



s— 1 

[n 



r(.- 



Fa 



2+j-s-' 



(8) 
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Obtaining the above form of the FC numbers, in which 
a ratio of products of the gamma functions of a shifted 
argument appears, is a key step of our reasoning. It 
allows us to represent the inverse Mellin transform of Eq. 
([8]) as a certain special function. To see this recall that 
the Meijer G-function of the ar gum ent z can be defined 
by the inverse Mellin transform pOj : 



G 



ini.n I „ 



(9) 



M- 



n™ir(ft+a)n;ur(i 



n 



ir(i-/3,-a)m^„+ir(a, + a)' 



This definition involves four lists of parameters, which 
can be represented by p complex numbers aj and other 
q complex numbers /3j. Integers numbers p and q can 
be equal zero and it s is assumed that < m < g and 
< n < p, so that possibly empty products in this form 
are taken to be equal to unity. A detailed description 
of the integration contours of the Mellin transform 
general properties of the Meijer functions and its special 
cases can be found in I40ll. 



P fx) 0.4- 




FIG. 1: (Color online) Marchenko-Pastur distribution Pi{x) 
compared with the Fuss-Catalan distribution P2{x). The sin- 
gularity at a; ^ is of the type Ps(x) ~ a;""'''''"''^'. 



of the argument z = xs^{s + 1)"*^'*+^). Looking at the 
range of the parameter j in ([8]) we see that the numbers 
of parameters of the Meijer G-function have to be set 
to n = 0, p = s and m = q = s. Hence this function 
involves 2s parameters, which read aj = (1 + j — s)/s 
and Pj = (j - 1 - s)/{s + 1) for j = 1, . . . , s. 

As there are only two products of gamma functions 
in ([8|), in contrast to four in (l9l), which is equivalent to 



P [x) 




10 12 14 16 
X 



FIG. 2: (Color online) Fuss-Catalan distributions Ps{x) plot- 
ted for parameter a equal to 3, 4, 5 and 6 are supported on the 
interval [0, (s-(- l)'''*'^/s'']. To show behaviour at the right side 
of the support the figure is plotted for a; > 5. 



setting n — Q and m = q, a further simplification of the 
above formula is possible. In this very case the Meijer 
G-function can be written as a combination of hyper- 
geometric functions of the same argument z - see Eq. 
(5.2.11) p. 146 of 



Direct comparison of expression ([8]) for the Fuss Cata- 
lan numbers and the Mellin transform as the Meijer G- 
function ([9|) allows us to represent the Fuss- Catalan dis- 
tribution Ps{x) by a Meijer G-function, 



Psix) = 



-3/2 



/2^ (S + 1)^+1/2 



G; 



/ I ai---ap 



(10) 



Let pi^g(^[{aj}^^J , denote the hyper- 

geometric function [42| of the type pFq with p 'upper' 
parameters aj and q 'lower' parameters bj of the argu- 
ment X. The symbol {ai}^^i represents the list of r ele- 
ments, ai, . . . Or. Then formula (fTO|) for the Fuss-Catalan 
distribution can be rewritten as 



Psix) 



E 

fc=i 



A 



k,s 



([{' 



ii,=iJ' Lr+s + i/j=i'i^+s + i/ 



j=k+i 



(s-t- l)''+i' 



(11) 



4 



where the coefRcients A^^^ read for fc = 1, 2, . . . , s 



,-3/2 



+ 1 /s-'/('*+^) X*^ 



27r 



s + 1 



n 



k-l 



n 



n 



S + l / 



(12) 



This formula, afong with Eg.ip^ befow, constitutes 
the key result of the present note. It gives an exact result 
for the Fuss-Catalan distribution Ps{x) for an arbitrary 
natural s. The FC distribution describes the density of 
squared singular values of a product of s independent 
square Ginibre matrices in the limit of large matrix size. 

The convergence conditions of the hypergeometric se- 
ries sFs-i immediately yield the support of Psix) which 
is equal to [0, (s + lY'^^/s'']. For small values of x the 
distribution behaves as x~^''^'^~^^\ It is contorting to see 
that in the simplest case s = 1 the above complicated 
form reduces indeed to the Marchenko-Pastur distribu- 
tion, 



„ , N 1 /r li r 1 1 \ \/l - x/4: , , 

Furthermore, the distribution P2(a^), shown in Fig. [TJ 
V3 



^^^")= 2;?^^^^(t"6'3J' 27" 



V3 



is equivalent to the form. 



6' 3^' ^3^' 27 
^3^ 27 



(14) 



^^3 [^(27 + 3V81-12x)^-6^] 
^2(2;) = —7;- — 1 , (15) 



127r 



;3 (27 + 3V81 - 12a 



valid for x E [0, 27/4] and obtained first in [3T| in context 
of construction of generalized coherent states from com- 
binatorial sequences. The distribution Psi^x), plotted in 
Fig. [21 is given by a sum of three terms. 



^3(2;) 



' 12' 4' 12J' i^' 256^ 



defined for n = 0, 1, ... . Here p and r are treated as 
integer parameters, p > 2, r — 1, 2, . . . . Setting r = 1 



and p = s + 1 we have ("^J 



(7ip+l\ _ 



fnp\ 



SO 



ri(p-l) + l V Ti 7 ' 

the numbers i?s+i^i(ri) are equal to FCs{n). Further 
relations involving the sequences Rp^ri^) are 



and 



Rp+i^p+i{n) = FCp{n + 1), 



Rp,pin) = Rp.iin + 1), 



(18) 



(19) 



which can be verified directly from their definitions ([T]) 
and (fT7| . 

The Raney lemma [43| implies that the number of the 
Rancy sequences of order p and length pn + 1 , for which 
all partial sums are positive, is given by the Fuss-Catalan 
numbers FCp-i(n) = i?p,i(n). Furthermore, as the num- 
ber of positive generalized Raney sequences is equal to 
Rp.rin) [111 we will refer to Rp,r{n) defined in p?| as 
Raney numbers. These numbers appear as coefficients 
in a generalized binomial series [22]. Some representa- 
tive examples of sequences Rp^rin), for n = 0,1,..., 7 
are quoted below together with their OEIS labels (3^ . 

RiA^) = 1,2,9,52,340,2394, 17710, 135720,. ..(A069271) 
i?5 2(n) = 1,2, 11, 80, 665, 5980, 56637, 556512, ...(v4118969), 

whereas two following sequences are not represented in 
OEIS: 

RiA'^) = 1,5,30,200,1425,10626,81900,647280,..., 
i?g3(n) = 1,3,21,190,1950,21576,250971,3025308,.... 



47ra;i/2 ^^^{Iq^ 2' 6^' ^'4^' 256 ^ 
^/2 i;, /r 5 3 13. .5 3. 27 \ 
"ei^^^n^I^'i'I^J' ^'2^' 256'^J- ^^^^ 

In Fig. 2 we present the Fuss-Catalan distributions Ps (x) 
for s = 3, 4, 5 and 6. 



III. RANEY DISTRIBUTIONS 

The Fuss-Catalan numbers FCs{n) defined in ([T]) can 
be considered as a special cases of a larger family of se- 
quences. 



Rp,r{n) := 



r f pn + r 



pn + r \ n 



(17) 



In a recent work Mlotkowski [20| has shown that the se- 
quence pT]) describes moments of a probability measure 
fip^r with a compact support contained in [0, 00), if point 
(r, p) which determines parameters of the Raney numbers 
belongs to the set S defined by inequalities p > and 
< r < p. Note that the point (1,1) implies a constant 
sequence of moments, i?i.i(n) = 1, which represents a 
singular, Dirac delta measure, — S{x — 1). 

In the case the measure fip^r is represented by a den- 
sity we will denote it by Wp^rix). Setting r to one 
one gets the Fuss-Catalan numbers, which implies that 
1^2,1(2:) represents the Marchenko-Pastur distribution, 
while Ws+i,i{x) reduces to the Fuss-Catalan probability 
density Ps{x). 

In general parameters p and r can be taken to be real 
and then the moments of the measure are expressed by 



5 



the gamma functions, 



T{np - 



np + r T(n + l)r(7ip + r — ?i + 1) ' 

(20) 

where the integration covers entire support of the mea- 
sure ^ip^r- For 1 < r < p the distribution Wp^ri^) is a 
positive function, see [20l | and Appendix B. 

The corresponding distribution Wp^r{x) can be written 
imphcitly by its 5'-transform, which allowed Mlotkowski 
to establish relations between various Raney distribu- 
tions [205, listed in Appendix C. For a precise definition 
of the 5-transform (or the free multiplicative transform) 
see Eq. 4.9 in [20]. In spite of these concrete results an 
explicit form of the Raney distribution Wp^r{x) has not 
appeared in the literature so far. 

Making use of the inverse Mellin transform and the 



Meijer function we can generalize results of the previous 
section and obtain explicit expressions for the Raney dis- 
tributions Wp^r{x), which correspond to integer values of 
the parameters p and r. 

Repeating steps analogous to Eqs.(|7l) - pO)) we can 
represent distribution Wp^r{x) in terms of the Meijer G- 
function. The explicit expression generalizing Eq. ljlOp 
reads 



Wr, 



P 



r-p-1/2 



QP.i 



Pi- 



(21) 



where the argument of the function is z = x{p— ly^^ /p^ , 
while its parameters are ai = 0, aj = (1 — P + 
for j = 2, . . . ,p and /3j = (r—p— l+j)lp for j = 1, . . . ,p. 
In analogy to pT|) this relation can be represented by the 
following sum consisting of of p terms: 



(p-l)P-l 



where pFp_i is the hypergeometric function and the numerical coefficients r2(p, r; j), for j — 1, 2, . . . ,p read 

^ 1 [nriir(^)l[nriir(^ 



'-^^^ V2^ (p - 1)'-P+'V2 V pP ) 



p—r+l-j \ Y[P J^f ' — P+ij _ r-p-l+j ' 
V p J lli3=2^V p-1 p 



(22) 



(23) 



Convergence properties of the hypergeometric func- 
tion imply that the Raney distribution Wp,r{x) for in- 
teger values of its parameters is supported in the interval 
[0,Kp_i], where is given in Eq. Formula 
implies that for p > r the distribution Wp^r{x) diplays a 
singularity for small x of the type x"'-^"'"-'/^. For p — r 
the 'diagonal' Raney distributions behave for small argu- 
ments as Wp^p{x) ^ x^Ip . 

It is helpful to add some clarifying remarks concerning 
the key formula . We draw attention to the fact that 
some simplifications will always appear for two reasons: 

a) firstly, one parameter form the 'upper' list of the 
parameters will be always equal to a parameter from the 
'lower' list. Consider, for instance, the case p = 4 and 
r = 2. Then the value of the first parameter in the 'upper' 
list becomes 1 + (j — 3)/4 and it cancels with the value of 
13 = 3 of the first sequence of the 'lower' parameters. One 
can demonstrate that a similar cancellation effect takes 
place for any pair (p, r). Therefore, the hypergeometric 
function pFp_i in (22) effectively reduces to p_iFp_2- 

b) secondly, we see from Eq. ([23l) that the coefficient 
^(Pi j) vanishes for j = p + 1 — r, due to the presence 
of the first gamma function in denominator. Therefore 
the sum in Eq. ([22l) involves (p — 1) terms with dif- 
ferent hypergeometric functions p_i-Fp_2- Then it can 
be explicitly verified that from the equality between the 



numbers i?s4-i,i(n) — FCs{n), the corresponding equal- 
ity between the probability distributions 

Ws+i,i{x) = Psix) (24) 

follows. In particular the following identity between the 
coefficients holds, r2(s-|-l, = Aj^s, whose demonstra- 
tion is rather tedious and will not be reproduced here. 

With the two provisos explained under items a) and 
b) above, Eg.ip^l '^iH used to obtain explicit forms 
of distributions Wp^rix) for small values of r and p. In 
particular, the case W2^2{x) reduces to the celebrated 
semi-circular law [5]: 

W,,2{x) = ^ iFo{[-l], [ ]; f ) = ^v/^(4^. (25) 

The above semicircle is centered at x = 2, while the 
Wigner semicircle centered at a; = is used in random 
matrix theory to describe the asymptotic level density of 
random hermitian matrices form Gaussian ensembles 0, 
The Raney distributions W2,rix) are plotted in Fig. [3] 
for r = 1 and r — 2. For comparison we have also plotted 
the function W2^3{x) furnishe d by Eq. ([^ . According 
to the results of Mlotkowski [20| and our Appendix B 
this function is not positive in its domain, so it does not 
represent a probability distribution. Here we obtain it 
explicitly. 
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/r, {x) 0. 

2, ' 





FIG. 3: (Color online) Raney distributions W^^rip^-) with 
values of the parameter r labeling each curve. For r = \ 
it reduces to Marchenko-Pastur distribution P\(x\ while a 
semicircle law is obtained for r = 2. For r = 3 the function 
represented by dashed line is not positive, which is implied 
by p < r. 



It is easy to observe that the above semicircle distribu- 
tion is related with the Marchenko-Pastur distributions 
P\{x) by a relation W2,2(a^) = xP\(x). This is a special 
case of a more general relation involving the "diagonal" 
Raney distributions with r — p and Fuss-Catalan distri- 
bution, 



FIG. 4: 



Wp,p{x) = X Wp^i{x) = xPp^iix). 



This result, established first in [20|, follows also naturally 
from Eqs.(lll|) and (HH). Thus for p = 3 one has 

Ws^iix) = P2ix), and W^^x) ^ ^ P2{x), (27) 

where the Fuss-Catalan distribution P2{x) is given in 
(ITil) . Due to an equivalent expression ([T5|) the distribu- 
tion Ws^slx) can also be expressed in terms of elementary 
functions. The intermediate Raney distribution, corre- 
sponding to r = 2, 



W3,2ix) 



V3 



Tj3 2-P'l(^ 



27rxi/3 --^VL 3' 6^' ^3^' 27^' 
V3a;i/3 /I 5 5 4 

^.^i([3-6]'y=^^) (28) 



I87 



27 



is shown in Fig. |4l In a close analogy to Eq. ([T5|) this 
distribution enjoys a similar representation in terms of 
elementary functions, 



W3,2ix) 



V3^ [(27 + 3V81 - 12a;) ' - 18^x1] 



367r 



x^ (27 + 3V81- 12a 



(29) 

Observe that the distribution 14^31(2;), W^,2{x) and 
^^3.3(2;) behave for small x as x"^/'^, and x^^^, 

respectively. 



(Color online) As in Fig|3] for Raney distributions 

31 



W3,r{x) supported in [0,6|]. The case W3^4,{x) (dashed line) 



is not a probability measure, 
value of r. 



Each curve is labelled by the 



Let US now discuss the case p = 4 illustrated in Fig. [S] 
Due to relations (l26l) one has 



14^4,1(2^)^^3(2;), and 1^4,4(2;) = 2; F3 (2;), (30) 

where the Fuss-Catalan distribution ^3(2;) is given in 
(|16p . Two intermediate Raney distributions have a sim- 
ilar form 



(26) ^4,2(2:) 



1 



3F2 



47ra;i/4 ^ 2 
^/2a;V4 



1111 3 

J' L4' 4J' 



1287r 



and 



^4,3(2:) 



V27ra;i/4 

_3V2x^ 

GAtt 
3.1/2 

3^2 



F2 



3F2 



6' 6' 2 

5 3- 
12' 12' 4- 
7 11 5 



27_ 
256 



([^ 



1 5 
2'4 
3 7 



12' 12' 4]' ^2' 4 



.) 

27 

256' 

i^i^)'<=") 



1 1 5 



-1 li 



27 



- 4' 12' 12^' ^4' 2^' 256 
1 7 11. r3 3. 27 



/ri (' 11-, r6 6-. 

^^U4'12'12^' U'2^' 



256 



327r 



-1 5 7. .5 7. 
■2' 6' 6^' ^'12^ 



27 ^ 
256' 



(32) 



The general formula ([22|) allows us to obtain an explicit 
form of the functions W4_5(x), which is not positive and 
it does not represent a probability distribution, see the 
dashed curve on Fig. 5. 

Figure [B] presents the Raney distributions Wp_r{x) in 
the "diagonal" case, r = p. As the distribution ^2,2(2;) 
represents Eq. ips]) . the diagonal Raney distributions, 
Wr^r{x), can thus by considered as a generalization of 
the semicircular law: they are defined for x £ {0,Kr~i) 
where Kg = {s + iy^^/s^, they are equal to zero at both 
ends of the domain and they are characterized by a single 
maximum. However, of r > 2 the functions Wr^rix) are 
not symmetric anymore, compare Fig. [Bl A plot of the 
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0.5- 




X 



-0.1 -I 

FIG. 5: (Color online) As in Fig. [3] for Raney distributions 
W4^r{x) . Dashed line corresponds to a quasi-measure ^4,5 (2;) , 
which is not positive. 




X 



FIG. 6: (Color online) 'Diagonal' Raney distributions 
Wr,r{x) form a generalization of the semicircle law 14^2,2 (a;). 
For small x they behave as x^^^ . Value of r labels each curve. 

parameter space {p, r) in which these distributions are 
marked together with the Fuss-Catalan distributions is 
presented in Fig. [T] 

IV. CONCLUDING REMARKS 

In this work we obtained an explicit form of the Fuss- 
Catalan distribution Ps{x). Obtained result is exact for 
an arbitrary s and it allows for a simple use of these prob- 
ability distributions. Results derived are relevant from 
the point of view of statistical physics as they describe 
asymptotic level density of normalized positive random 
matrix of the product form X — (Gi • • • Gs){Gi ■ ■ ■ GsY , 
where Gi, . . . ,Gs denote s independent random matri- 
ces from the complex Ginibre ensemble. The variable 
X = NX denotes the rescaled eigenvalue A of random 
matrix X of size TV. 

It should be emphasized here that the Marchenko- 



r 




1 2 3 4 5 6 



FIG. 7: Raney distributions in the parameter space {p,r): 
Marchenko-Pastur distribution Pi is represented by the point 
(2, 1), while Fuss-Catalan distribution Ps is represented by (*) 
at {s + 1, 1). The semicircle (sc) law corresponds to (2, 2) and 
the 'dia gon al' Raney distributions (-I-) to (p,p). As shown 
in Ref. [2^ any point in the shaded set E corresponds to a 
probability measure with a compact support in [0, 00). 

Pastur and Fuss-Catalan distributions describe the level 
density of Wishart like random matrices in the limiting 
case N ^ 00 only. In practice, for any fixed A'^ the finite 
size effects occur. As discussed by Blaizot and Nowak 
[3 SI] the finite N effects are related with the diffraction 
phenomena, while the large N limit of the random ma- 
trix theory may be compared with the geometric limit of 
the wave optics or the semiclassical limit of the quantum 
theory. An explicit description of finite N corrections to 
the spectral density of Wishart matrices obtained from 
products of Ginibre matrices is provided by Burda et al. 

A simple argument put forward in [18| (see Appendix 
B therein) implies that the same FC distributions de- 
scribe also the spectral density of a product of s random 
matrices taken from the reai Ginibre ensemble [47j . The 
case of a product of two real matrices, recently studied 
in context of quantum chromodynamics 48] , is also im- 
portant in applications in econophysics, where one uses 
a product of two real correlation matrices [T^ . 

The Catalan numbers are a special case of a more 
general, one parameter family of Fuss-Catalan numbers, 
which from a subset of two-parameter family of Raney 
numbers. In the same way, the Marchenko-Pastur dis- 
tribution -Pi (a;) is a special case of the Fuss-Catalan dis- 
tributions, which in turn belong to the two parameter 
family of Raney distributions Wp^rix). This wide class 
of probability distribution includes e.g. the Dirac delta, 
6{x — 1) and the semicircle law, ^2,2(2;). Applying the 
inverse Mellin transform for integer parameter values of 
the parameters we found an explicit exact representation 
of the Raney distributions in terms of the hypergeometric 
function. For any r — 1,2, ... ,p the Raney distribution 
Wp,r{x) is supported in the interval [0,p^/{p— 1)^"^]. 
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For small x the distribution behaves as 



if r < p, 
if r = p. 



(33) 



The Raney numbers (|17p imply that the mean value of 
the Raney distribution Wp^r{x) is equal to r, while the 
second moment reads r{2p + r — l)/2. 

Let us conclude the paper with the following remark. 
The Fuss-Catalan distribution describe statistical prop- 
erties of singular values of products of random matri- 
ces of the Ginibre ensemble. It is then natural to ask, 
whether there exist any ensembles of random matrices, 
such that their squared singular values can be described 
by the Raney distributions. 
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The series (IA2[) . if convergent, defines a generalized 
hypergeometric function, which may then be defined over 
a wider domain of the argument by analytic continuation. 

Note that for s = 2 the expression (fTT|) involves the 
Gauss (ordinary) hypergeometric function 

2Fi([a,6],[c];z) = f;^%i^^, (A3) 



that includes many other special functions as special or 
limiting cases. 



Appendix B: Positivity of the distributions 

Psix) and Wp,r{x) 

For completeness we prove in this appendix that for 
any integer s the distribution Ps{x) is positive for x G 
{0,K,), where Ks = {s + 1)*+Vs". Eq. © implies that 
Ps{x) is given as the inverse Mellin transform of Eq. ([8]). 
We are going to use the convolution property for two 
Mellin transforms, M\f (x); a] — f*{(j) and M[g{x); a] — 
g*{(j), which reads [36|, |37[ 

M-'[na)g*{a);x]^ f{^)g{t)j 



dt 



(Bl) 



Appendix A: Some information on special functions 

Special functions of interest for this note are related 
to the Euler gamma function, which admits a following 

integral representation, T{z) = J V'^^ e^'^ dt. Integrat- 



ing by parts we see that V{z -1-1) — z ■ T{z). For an 
integer argument the Euler function is given by factorial, 
T{n + 1) = n\ . The Gauss-Legendre formula allows one 
to compute the Euler gamma function of a multiple of 
argument [i^ . 



fe-i 



r(fcz) = (27r)(i-'=)/2 -Q ^ U (A^) 

i=0 



for fc = 1, 2, 3, . . . and z 7^ 0, -1, -2, . . . . 

The generalized hypergeometric series is a series in 
which the ratio of successive coefficients indexed by n 
is a rational function of n. It can be defined as 

(A2) 

where we use the Pochhamnier symbol, defined by (a)„ = 
a{a + l)(a -I- 2). ..(a + n — 1) and (a)o = 1- 



If X > and both functions f{x) and g{x) are positive 
then their Mellin convolution defined by the integrals 
(|B1[) conserves positivity. 

Consider now, for a given j = 0, 1, . . . , (s — 1), the in- 
dividual term in the product in Eq.®. Its Mellin trans- 
form will satisfy, due to formula 8.4.2.3, p. 631 in [i^l, the 
following equality 



r T a 



s+1) 



2+.;-s 



'(l-x) 



' 2{s+l)+js-s^ - 
s+1 



(B2) 

which for all j = 0, 1, . . . , (s — 1), is a positive function 
for X G (0, 1). Then Eq. ^ can be viewed as a (s — 1)- 
fold convolution of positive functions, which by (IB1[) is 
itself positive for x G (0, Ks). The upper edge Kg of the 
support can be read off from the prefactor in Eq.®. 

In a similar way one can prove positivity of Raney 
distribution Wp^rix) for natural values of the parameters, 
provided r < p. Here is the streamlined version of the 
proof. We use the aforementioned formula 8.4.2.3 of [i^l 
for b > a, 



Tia- 



r(a + h) 



\r{b-a)\ 



\b-a-l 



(B3) 



which describes a positive function for < a; < 1. We 
quote now the full version of the analog of Eq.® for 
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Rp.r{<7), with (7 = n + 1 and p = 2,3, 



P 



,r-p-l/2 



r P'' - 
r(a + ^:^) r(a + ^:^£±i 



r(a) 



.(B4) 



Consider first the case 1 < r < p. The weight func- 
tion Wp^r{x) is the inverse Mehin transform, Wp^r{x) = 
M~'^[Rp,ricr);x], and, from ((Bit via Eq. pTI) . it is the 
p-fold Melhn convolution of r[(T + (r — p)/p]/T{a) and 
of p — 1 factors in the product in Ea. (jB4l) . For each of 
these individual ratios of gamma function relation (jB3p 
holds. Therefore, with Ea. (|B4p . the positivity of the dis- 
tribution Wp^r{x) follows. 

In the second case, r ^ p, only p—1 factors in Ea. (jB4[) 
intervene in the convolution. Hence Wp^p{x) is also posi- 
tive, li r > p the first ratio in Eq. (jB4l) destroys the pos- 
itivity, so no further considerations are needed to show 
that in this case the function Wp^r (x) is not a probability 
distribution. 



Appendix C: Relations between Raney distributions 



The Raney distribution Wp^rix) may be implicitly 
written by its S'-transform, which allowed Mlotkowski 



[20| to establish a relation between distributions with 
various values of their parameters in terms of the free 
multiplicative convolution denoted by (compare Eq. 
4.10 in 201), 



(CI) 



This result, valid for p,q > 0, can be generalized for an 
arbitrary positive r. 



p-\-rq,r ■ 



(C2) 



Moreover, there exists yet another relation. 



[W: 



H-P.iJ 



sp,l 



(C3) 



which holds for s > and is equivalent to the multiplica- 
tive convolution property for the Fuss-Catalan distribu- 
tion. 
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